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We study Kondo physics in the algebraic spin liquid, recently proposed to describe 
ZnCu 3 (OH) 6 Ch [Phys. Rev. Lett. 98, 117205 (2007)]. Although spin dynamics of the alge- 
braic spin liquid is described by massless Dirac fermions, this problem differs from the Pseudogap 
Kondo model, because the bulk physics in the algebraic spin liquid is governed by an interacting 
fixed point where well-defined quasiparticle excitations are not allowed. Considering an effective 

£ — ■ bulk model characterized by an anomalous critical exponent, we derive an effective impurity action 

in the slave-boson context. Performing the large- JV CT analysis with a spin index N a , we find an impu- 
rity quantum phase transition from a decoupled local-moment state to a Kondo-screened phase. We 
evaluate the impurity spin susceptibility and specific heat coefficient at zero temperature, and find 
— . ■ that such responses follow power-law dependencies due to the anomalous exponent of the algebraic 

spin liquid. Our main finding is that the Wilson's ratio for the magnetic impurity depends strongly 

jy^ . on the critical exponent in the zero temperature limit. We propose that the Wilson's ratio for the 

magnetic impurity may be one possible probe to reveal criticality of the bulk system. 

m ■ 

i—l | PACS numbers: 71.10.-w, 71.10.Hf, 71.27. +a 

\ I. INTRODUCTION 

Recent experiments have claimed the emergence of spin liquid (SL) phases in materials of geometrically frustrated 
lattices such as Cs2CmC74,[I[ k— (ET)2Cu2(CN)3,^ and ZnCu^{OH)QCl2^, where no symmetries associated with 
spin rotations (magnetic ordering) and lattice translations (valance bond ordering) are broken at low temperatures 
while charge fluctuations are frozen due to strong electron-electron interactions (Mott insulator). An important issue 
is the nature of such SL phases. Although spin susceptibility, specific heat, and thermal transport measurements can 
<— | ■ determine possible spin liquids, there still remains uncertainty. 

Consider CS2C11CI4 with an anisotropic triangular lattice. [1] Although this material exhibits magnetic long-range 
O • spiral ordering below T = 0.62K with an incommensurate wave vector, the spin-fluctuation spectrum in inelastic 
neutron scattering experiments has shown large high-energy continuum beyond the spin- wave description. In addition, 
£S) . this continuum spectrum survives above the Neel temperature. More detailed analysis revealed that the continuum 
' follows Imx(uj) ~ uj~ v with an anomalous exponent 77, suggesting the presence of dcconfined critical spinons. Such 
spin-fluctuation measurements suggest several candidates of SL scenarios, for example, decoupled one-dimensional 
"^t" ■ chains, Q proximate gapped SLs,Q algebraic spin liquid (ASL),[|| algebraic vortex liquid, [?J and so on.Q 

\ Recently, Florens et al. studied the role of magnetic impurities in both the Z2 SL phase and the 0(4) QCP 
* . ■ separating the spiral magnetic order from the Z2 SL.Q Although impurity moments coupled to spin-1 bosons (spin 
singlet-triplet excitations) in conventional paramagnets are only partially screened even at the bulk 0(3) QCP,[lO( 
they have shown that the presence of deconfined bosonic spinons can display a bosonic version of the Kondo effect. 
Furthermore, they found a weak-coupling impurity quantum phase transition (I-QPT) from a local-moment state 
to a fully-screened phase. This study implies that the magnetic impurity can be utilized as a probe for elementary 
excitations, thus identifying the nature of SLs. 

In this paper we investi gate the Kondo effect in the ASL, recently proposed to be realized in the Kagome antiferro- 
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magnet ZnCu 3 (OH) e Cl2, 11] where no magnetic order is observed down to very low temperature 50mK compared to 
the Curie- Weiss temperature (> 200K), and there is no sign of spin gap in dynamical neutron scattering. [i| However, 
it is not perfectly clear whether all experiments are consistent with the ASL conjecture. The ASL picture is not 
consistent with the temperature-linear specific heat below 0.5K and saturation of the spin susceptibility to a finite 
value below 0.3if ,[ll| because these measurements indicate the existence of a finite density of states at the Fermi 
energy. This discrepancy may result from the presence of disorder in real materials. To examine the role of magnetic 
impurities in the ASL can be an important test in revealing the genuine nature of the SL phase of this compound. 

The ASL can be found from the fermion representation of the Heisenberg model via the flux mean-field ansats.[l2j 
This reminds us of the previous study of the Kondo effect in the flux phase by Cassanello and Fradkin. 1131 More 
generally, one may regard the present impurity problem as the class of the Pseudogap Kondo model.(T3,[l3.ll5j where 
the fermion density of states vanishes as p(e) ~ |e| r near the Fermi energy. The case of r = corresponds to Fermi 
liquid while the r = 1 case coincides with Dirac fermions arising from the flux phase or d — wave superconductor. 
In contrast with the Kondo effect of the Fermi liquid, the Pseudogap Kondo model has shown that Kondo screening 
of the magnetic impurity can appear beyond some critical value of the Kondo coupling constant. Thus, the I-QPT 
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from a local-moment state to a Kondo-screened phase was found in this model. Furthermore, the exponent r in the 
density of states was shown to play the role of an effective dimension in the problem. The r = 1 case was found to be 
its upper critical dimension, thus exhibiting logarithmic corrections to scaling while the case of r = lies in its lower 
critical one. 

However, there is an important difference between the Pseudogap Kondo problem and ASL Kondo physics. The bulk 
physics in the Pseudogap Kondo problem is governed by a noninteracting (Gaussian) fixed point, thus allowing well- 
defined electron-like quasiparticle excitations. On the other hand, the ASL physics is determined by an interacting 
fixed point (the conformal invariant fixed point of QED.^),[l6| where well-defined spinon quasiparticle excitations 
corresponding to electrons do not exist. The absence of quasiparticle excitations prohibits us from applying the 
conventional picture of the Pseudogap Kondo physics to the ASL Kondo problem. In this respect the Kondo effect 
at such an interacting fixed point is an interesting problem. 

The main difficulty is how to introduce the absence of well-defined spinon excitations in the ASL Kondo problem. 
Long-range gauge interactions would result in the anomalous critical exponent rj^p in the single spinon propagator, 
destroying the quasiparticle pole in the Green's function. Unfortunately, such critical physics can be found within 
the summation of infinite diagrams of gauge interactions, and this procedure prohibits us from analyzing the ASL 
Kondo problem in a simple mean-field way such as the large- N a approximation with a spin index N a , well utilized 
in the Pseudogap Kondo problem. [l3j Considering the mathematical derivation in the large- N a context, the main 
problem is how to derive an effective impurity action from the ASL-Kondo Lagrangian through integrating out bulk 
degrees of freedom, critical spinon and gauge fluctuations coupled to the magnetic impurity. More precisely speaking, 
a bulk-spinon propagator appears to govern the impurity dynamics in the effective impurity action, thus how to write 
its accurate form is an important problem since the presence of gauge interactions makes such a task nontrivial. 

In the present paper we assume the expression of the spinon Green's function as an ansatz, introducing an anomalous 
critical exponent rj^. In the text we discuss the validity of this ansatz in great detail. This effective representation 
allows us to analyze the ASL Kondo problem in the large- N a context. Performing the slave-boson saddle-point 
analysis for the effective impurity action, we find an I-QPT from a decoupled local-moment state to a Kondo-screened 
phase. We evaluate the impurity spin susceptibility and specific heat coefficient at zero temperature, and find that 
such responses follow power-law dependencies due to the ASL anomalous exponent. The main finding of the present 
study is that the Wilson's ratio for the magnetic impurity depends strongly on the ASL critical exponent in the zero 
temperature limit. We propose that the Wilson's ratio for the magnetic impurity be a probe to reveal criticality of 
the bulk system. 

II. REVIEW OF THE ALGEBRAIC SPIN LIQUID AND ITS KONDO PROBLEM 

For completeness of this paper, it is necessary to review how the effective Lagrangian so called QED3 describing the 
ASL is derived from a microscopic model such as the antiferromagnetic Heisenberg model, H = Y)^ JijSi-Sj with Jy > 

0. Inserting the fermion representation of spin Si — \ fiercer' ha 1 m t° the Heisenberg model, and performing 
the Hubbard-Stratonovich transformation for an exchange channel, we find an effective one-body Hamiltonian for 
fermionic spinons (f ia ) coupled to a hopping parameter {xij), H e f f = -J2ija JijfiaXijfj(r+12ij Jij\Xij\ 2 - Notice that 
the hopping parameter Xij is a complex number defined on links ij. Thus, it can be decomposed into Xij — \Xij\ e j > 
where \xij\ and 9^ are the amplitude and phase of the hopping parameter, respectively. Inserting this representation 
of Xij hito the effective Hamiltonian, we obtain H e ff = —^2 i j rT Jij\Xij\fia eieii fjiy^ where the constant contribution 
for the ground state energy is omitted. Then, we can see that this effective Hamiltonian has an internal U(l) 
gauge symmetry, H' ef f [f^,9' i:j ] = H e ff[f i<T ,9ij] under the following U(l) phase transformation, f' la = e**« f lr7 and 
= @ij + 0j ~ 4>j- This implies that the phase 0y of the hopping parameter plays the same role as the U(l) gauge 
field ay. 

One can perform a saddle-point analysis of the effective Hamiltonian to find its stable mean-field phases in various 
lattices such as square, [l2| triangular, @ Kagome,[ll], [TtJ and etc. In the present paper we consider the square lattice 
for simplicity, where the antiferromagnetic long-range order can be suppressed via next-nearest-neighbor or ring- 
exchange interactions causing frustration. It is not so difficult to extend the mean-field analysis on the square lattice 
into that on the Kagome lattice, proposed to show the SL physics of ZnCui(OH)s\CU .[If| 

It has been shown that one possible stable mean field phase is a 7r-flux state, where a spinon gains the phase of it when 
it turns around one plaquette. The amplitude of the hopping parameter is frozen to be \xij\ = I (/jo- /«£>") I = Xo 
in the low energy limit. Then, one finds the low-energy effective Lagrangian in terms of massless Dirac fermions 
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interacting via compact U(l) gauge fieldspi 
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o-=T,|n=l 

Here, -071(7 is the two-component massless Dirac spinon, where n = 1,2 represent the nodal points of (ir/2,ir/2) and 
(7r/2, — 7r/2), and cr =|, L SU(2) spin. They are expressed as V'Ict = I ^ le<T ) an d V^er = ( '{ 2ocr ), respectively. 

V Jloe J \ J2ea J 

In the spinon field / n ; CT n = 1,2 represent the nodal points, I — e,o, even and odd sites, and a =t, I, its spin, 
respectively. The Dirac matrices are given by the Pauli matrices — (03, 02, cri), satisfying the Clifford algebra 
[7M;7f]+ = 2(5^^. is the U(l) gauge field whose kinetic energy results from particle-hole excitations of high energy 
spinons. e is an effective internal charge, not a real electric charge. 

It has been argued that QED3 has an infrared stable fixed point showing the conformal symmetry in the large- N a 
limit (cr = 1, iV CT ).[16| This conformal invariant fixed point is identified with the ASL, displaying algebraically 
decaying correlation functions with anomalous critical exponents. To confirm the ASL as a genuine stable phase a 
cautious person may ask the stability of such an interacting fixed point against perturbations. Four-fermion interaction 
terms are irrelevant at this fixed point owing to their high scaling dimensions. In addition, chiral symmetry breaking 
due to noncompact gauge fluctuations has been shown not to occur in the Schwinger-Dyson-equation analysis when 
the flavor number of massless Dirac fermions is sufficiently large. [l9[ Furthermore, it has been argued that confinement 
as an instanton effect arising from compact gauge fluctuations does not seem to appear in the large- N a limit because 
the scaling dimension of the monopolc insertion operator is proportional to the flavor number N a , thus expected to 
be irrelevant in the large- N a ASL. [la] 

Criticality of the ASL is characterized by critical exponents of correlation functions. The single particle propagator 
GASL{k) — (t/w(&)Vw(&)) can be expressed as 

G *SL(k) « -i^, (2) 

where r]^ is an anomalous critical exponent. One can find such an anomalous dimension in the large- N a analysis. [20i| 
However, the critical exponent obtained in this way is difficult to have a definite physical meaning because it is not 
gauge invariant. In this respect the critical exponent rj^ should be evaluated in a gauge invariant way. The following 
gauge invariant Green's function can be considered, Gasl(x) = (T T [ip n a(x)e l -fo ^f^^VwdO)]). Unfortunately, it 
is not easy to calculate the critical exponent with such a gauge invariant expression. Its precise value is far from 
consensus and still under current debates. The crucial point is the sign of the exponent fty while its absolute value is 
given by 1 77^, I ~ N^ 1 in the 1/N„ approximation. (20j . Most evaluations [2ll I22I I23I |24| suggest its negative sign, 77^, < 0. 
However, as argued in Ref. [25|, its negative sign seems to be unphysical in the sense that the spinon propagator 
becomes more "coherent" at long distances than the propagator of the free Dirac theory. This result is in contrast 
with the usual role of interactions, making elementary excitations less coherent. This is indeed true in such critical 
field theories with local repulsive interactions, for example, the iV-vector model, where positive critical exponents are 
well known. [2^| If the critical exponent is positive, long-range gauge interactions destabilize the quasiparticle pole. 
The quasiparticle weight Z(p) ~ p n ^ with momentum p vanishes in the long- wave length and low-energy limits. In 
the present paper we do not determine its sign. Instead, we regard the exponent as a phenomenological parameter. 
Thus, we consider both cases of 77^, < and 77^ > 0. Furthermore, we assume that the renormalized spinon propagator 
[Eq. (2)] is obtained in a gauge invariant way. [2ll [22I [23l. [24l. [25| and the critical exponent -q^ is also gauge invariant. 

Another important character of the ASL is that the conformally invariant fixed point has an enlarged global sym- 
metry beyond the original lattice model, here the Heisenberg Hamiltonian. Such an emergent symmetry corresponds 



to Sp(4) in the case of SU(2) gauge interactions [271] and SU(4) in the case of U(l) ones[28|. This enlarged symmetry 
gives rise to an important effect on correlation functions, that is, resulting in the same behaviors between different 
correlation functions when the operators in the correlators are related with symmetry transformations. For example, 
staggered spin correlations have the same functional dependency (power-law decay) as the valance bond fluctuations 
since they are symmetry-equivalent. An interesting point is that such correlations are most susceptible in the ASL.(28| 
This implies that the ASL resides near the antiferromagnetic and valance bond solid phases. Actually, Tanaka and 
Hu have derived an effective Wess-Zumino-Witten (WZW) Lagran gian from the ASL, describing competition between 
antiferromagnetic spin correlations and valance bond fluctuations. [29( 

To study the role of magnetic impurities in the ASL bulk, we consider the Kondo coupling term, Hk = % Sq'S, 
where S q is a spin-fluctuation operator of bulk spinons with momentum q and s represents an impurity spin. The 
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bulk-spin operator has two contributions in the continuum, 

S(q) ~ S u (q) + S s (q) = E E ~ <?)7o^y-?/w(fc) + i> n *(k - g)^-Vw(&)> (3) 

k n<T<r f k naa' 

where S u (q) represents the uniform component and S s (q) denotes the staggered one.[l8j Then, the ASL Kondo 
problem is described by the following action 

2 

S = dr{ d 2 r( J2 E ^7„(^ - ia^na + X\d x a| 2 ) + ^ ^(S u (q) + S s (q)) ■ s). (4) 
o-=T4™ =1 ' i 

The next work is to obtain an effective impurity action, integrating out bulk degrees of freedom, spinon and gauge 
excitations coupled to the magnetic impurity. One can write down its schematic expression in the following way 

J q q 

where (^^(g, t)S^ u )(— <7, r')) is the renormalized correlation function of staggered (uniform) spin fluctuations, and 
• • • are higher moment contributions. As clearly shown in this expression, dynamics of impurity spin fluctuations 
is governed by spin correlations of the bulk at the impurity site. An important point is that only staggered spin 
correlations exhibit an anomalous scaling behavior with a nontrivial critical exponent. 1 30, 31] Uniform spin correlations 
have no anomalous scaling dimension since they correspond to conserved currents. [18l l30j Correlations of conserved 
currents do not have any anomalous scaling dimensions. This means that the contribution of uniform spin fluctuations 
is basically the same as the Kondo effect of the Pseudogap Kondo model while that of staggered spin excitations will 
give rise to new effects on the Pseudogap Kondo physics. Furthermore, staggered spin fluctuations are most singular 
in the large- N a ASL, [5^] thus expected to contribute to the Kondo effect dominantly. In this respect we take into 
account staggered spin fluctuations only, which is an important assumption in the present paper. 



III. KONDO PHYSICS IN THE ALGEBRAIC SPIN LIQUID: LARGE- iVa ANALYSIS 

Our objective is to construct a mean-field theory for the present Kondo problem. Using the slave-boson rep- 
resentation, the impurity spin is expressed as s = | ^ , Xa^crcr 1 Xa' > an d such fermions satisfy the constraint 



J2a=i XaXa = Q x with Q x = 2s, where s is spin. Inserting this expression into Eq. (4) with Eq. (3), the Kondo 
coupling term becomes 

H K = E E E $nv(k-q)XvXi'10lpn<T>(q)- ^jy^EE E ^( fc ~ q)Xaxi^na'{q) (6) 

a qk n— 1 (jo-' — 1 G qh n—laa'—l 



in the large- N a treatment, where the first and second terms are associated with uniform and staggered spin-fluctuation 
contributions, respectively. Since staggered spin fluctuations will give main contributions to the ASL Kondo effect, 
effects of uniform spin fluctuations arc neglected in the following. 

Performing the Hubbard-Stratonovich transformation for the Kondo-exchange channel, we find an effective ASL- 
Kondo action as our starting point 

Seff = J d 3 x[{j2J2^n^(d^ -ia^na + x a| 2 ) -^^K^tfOl+tClM^:) 

(T—1 n—1 a—I n—1 

, N, iV„ N n 

+ j dr[Y^ xUdr ~ K) Xa + xUa - Q X ) + E b n^ b n] ■ ( 7 ) 



a— 1 n—1 



The first part represents the ASL bulk. The second part arises from the Hubbard-Stratonovich decoupling of the 
Kondo interaction term, where b s n is a two-component hybridization order parameter associated with staggered bulk- 
spin fluctuations. Such a hybridization order parameter is determined self-consistently in the saddle-point analysis 
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2J K ' u ™ V (2T)Vi 
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The third part describes impurity-spinon dynamics, where h a = aH is an external magnetic field and A is a Lagrange 
multiplier field to impose the pseudo-fcrmion constraint. 

When the bulk system is in the non-interacting fixed point corresponding to the absence of gauge interactions, the 
effective Kondo model becomes the multi-channel Pseudogap Kondo model, where the channels come from Dirac nodes 
n = 1, ...,N n . This model was argued to show an I-QPT from a decoupled local-moment state to an over-screened 
phase in the large- N a approximation although this analysis does not capture the over-screened Kondo physics quite 
well.[l3| On the other hand, the present bulk system lies at the interacting fixed point characterized by the anomalous 
critical exponent r)^, where quasiparticle excitations do not exist. In this case it is not clear whether the conventional 
Kondo screening picture is applicable. 

Integrating out bulk-spinon and gauge excitations, we obtain an effective impurity action in energy-momentum 
space 

K7f = / 3Exl(*-i+e x k,-EE^( 44<vw(fc)^(fc)))^x^ 

9J -T, b n^b s n -e x Q x }, (9) 
2Jk t^i J 

where iX is replaced with e x to clarify its physical meaning. The main question in this impurity action is how to evaluate 
the spinon Green's function. As discussed intensively in the previous section, the single particle propagator has an 
anomalous scaling exponent, given by (Vw(&)Vw(&)) = — ijfik fJ ,/\k\ 2 ~ rli ' . This expression seems to be consistent 
with Eq. (5) if " sub- leading" uniform spin-correlation contributions are not taken into account. This is because 
the critical exponent of the staggered spin-spin correlation function is found to be twice the exponent of the single 
particle propagator, i.e., 2r]^ in the case of rj^, < 0.[13,[3l[ Such correspondence occurs when both critical exponents 
are calculated in a gauge invariant way. This correspondence was also pointed out in Ref. HH 

Such spinon excitations with an anomalous scaling exponent result in anomalous energy-dependent (nonlocal in 

time) interactions for impurity fermions, as reflected in the kernel of J i^pM^ = ij^F^ko)- The vector function 
Ffj.(ko) is obtained to be 

f d 2 k k, (fcg + A^-ifcpi^ 

F ^ 0) ~ J (2n) 2 (Jfeg+ifc2)l-W2 " fc ° 2^ (10) 

thus ij^F^ko) = ijokoF(ko) with F(ko) = [(fc 2 + A 2 )* 1 '*/' 2 — |fco|' 7 *]/27i7fy,, where A is a momentum cutoff. 

Inserting Eq. (10) into Eq. (9) and integrating over impurity fermions in Eq. (9), we find the following expression 
for the impurity free energy 



F — 

L imp — 



^ E ln [(* fc ° -K + e x ) + ik F(k ) ]T btfb>\ + ^ E b '^oK - e x Q x . (11) 



n— 1 n— 1 



Expressing the hybridization order parameter as a two-component spinor bfj = (&„ + 6„_), one can find b s n _ — in the 
saddle-point analysis. Representing the above impurity free energy with b s n+ = 26, we obtain 



T? - r dk ° 

mp ~ 2 J_„ 2n 



In 



InUko -H + e x + l k ^\b\ 2 {(k 2 + K 2 )^' 2 - \k Q \^}) 
2N n 



(ik + H + e x + iko '^\b\ 2 {(k 2 + A 2 )W2 _ | feo | W })] + ^^| 6 | 2 _ e%Q (12) 



Minimizing the impurity free energy with respect to b and e x , we find the saddle-point equations giving the self- 
consistency 



1 r dk k 2 ^m + a 2 )"- /2 - « [1 + t&ww + a 2 )W 2 - w }]. 

(13) 



:j o ' / j L -Krj-ip I I \ v 

-r ii ; jvq jj v 

VJk Jo ~2^ k 2 [l + ^\b\ 2 {(k 2 + A2)^/2 _ fc^'}] 2 + e\ > ~ 

! X _ 1 dk 2e x 



X " Jo 27r fco 2 [l + ^l&l 2 {(fco 2 + A 2 )^/2-fc^}]- 



Since the impurity free energy is momentum-cutoff-dependent, it is necessary to make it cutoff-independent, taking 
appropriate scaling transformations for all variables. Considering the scaling dimension of ij) ncr given by dim[^„ CT ] = 
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FIG. 1: Renormalized critical Kondo coupling constant as a function of the ASL exponent. 



1 + r}ifi/2, one can find dim[6] = —rj^,/2 and dim[Jjc] = —1 — T]^, where dim[0] represents the scaling dimension of an 
operator O. Then, the scale-free impurity free energy is obtained to be 



F, 



imp 



imp 



A 



- — — I il.r 

4tt 



2N 

ln( ix — h + e r + ix \b r \ \{x 



W 2 



+ ln(ix + h + e r + ix— \b r \ 2 {{x 2 + l)^/ 2 - \x\n\ 



— IM ~ < 



(14) 



where such rescaled variables are given by 

b 



b r = 



Jk 



ko , H 



Notice that these scaled variables are dimensionless. Accordingly, the self-consistent saddle-point equations read 



1 

Ox 



1 

77 .'0 



, sl{(x 2 + l)^ /2 - x^}[l + f^\b r \ 2 {(x 2 + 1)W 2 - X^}] 
dx— — 



dx- 



c 2 [l + ^\b r \ 2 {(x 2 + 1)W 2 - a?TJ*}]' 



T Jo X 2 [1 + ^\W\ 2 {(X 2 + 1)W2 + , 



(15) 



The QCP of the I-QPT can be found with b r — > and e r — > in the particle-hole symmetric case, Q x /N a = 1/2. 
Then, the critical renormalized Kondo coupling constant is obtained from Eq. (15), 



1 



1 



J, 



rc V4> JO 



dx{(x 2 + l)w/ 2 -x**} 



2tt 3 /2 w r(-at) ' 



(16) 



as far as the ASL exponent lies in — 1 < < 1. In the large N a limit the ASL exponent may also satisfy this 
condition as discussed in Section II. In addition, this critical value is continuously defined in the limit of r/^, — > ±0, 
where the impurity critical point [32J is given by 



J' r< 



2tt 2 



dx ln( 1 



1 



1 

2^' 



(17) 



consistent with the previous study, [la] 

It is interesting to notice that the I-QPT occurs as long as the ASL exponent |r^| < 1. Remember that in the 
regime of < r/^ < 1 critical spinon excitations are less coherent than those in the Pseudogap Kondo model = 0) 
while in the regime of —1 < < such spinon excitations become more coherent than quasiparticle excitations in 
the Fermi liquid with pseudogap. To screen the magnetic impurity, stronger Kondo couplings would be required when 
quasiparticle excitations are less coherent. Actually, we find such an asymmetric behavior for the ASL exponent in 
Fig. [TJ obtained from Eq. (16). 
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(a) (b) 

FIG. 2: (Color online) Impurity chemical potential, e r , and hybridization amplitude, b r , (log- log plot) as a function of the 
renormalized Kondo coupling, J r , for various ASL exponents. Here e r , b r , and J r are dimensionless rescaled variables. 



It might seem mysterious that the critical Kondo coupling vanishes as 77^ — ► ±1. As the ASL exponent approaches 
1, critical spinon excitations are not only less coherent but also localized. Considering the spinon propagator Eq. (2), 
77,/, = 1 makes it energy- momentum- independent. Such localized spinons are expected to form a Kondo singlet with 
an impurity spin immediately. When the ASL exponent goes to — 1, it is important that the bare scaling dimension 
of the Kondo coupling (dim[JA'] = —1 — Tty) vanishes, implying that Kondo interactions are marginal perturbations 
similar to the conventional Kondo effect in the Fermi liquid. In this respect the critical Kondo coupling would go to 
zero as — ► — 1. 

Solving Eq. (15) numerically, one can find the hybridization amplitude \b r \ 2 as a function of the Kondo coupling 
J r . We show the I-QPT in Fig. 2, where both b r and e r vanish as J r — > J rc . It is important to notice that the cc-axis 
is J r — J rc instead of J r . This means that the impurity QCP matches the origin of the £-axis. The absolute value 
of the impurity chemical potential increases rapidly as the ASL exponent increases from j]^ = —0.2 to rj^p = 0.2 [Fig. 
Ufa)]. Accordingly, the increasing ratio of the hybridization order parameter is largest for 77^ = 0.2 and smallest for 
??0 = —0.2. This may be associated with localization tendency emerging from a positive exponent. A further analysis 
finds a scaling behavior of the hybridization amplitude not only near the impurity QCP, but also rather away from 
the QCP, i.e., in the Kondo-screened phase. Such a scaling behavior even in the Kondo phase seems to arise from the 
criticality of the bulk system. From the log-log plot of Fig. E[b), we nn d the scaling relation 

I6 r | 2 ~ (j r - j rc y^ (is) 

with f(r)i/)) ~ 3 + 2?7^, confirming that the slope of the positive ASL exponent is larger than that of the negative one. 
The I-QPT can be also found in the impurity-spin susceptibility, 

gWfe) = N a r°° dx ^-x a [i + ^IM 2 {(^ + i)^ 2 -^}] 2 

9h 2 TT J X ( x 2[ 1 + M^| 6r |2|( a; 2 + 1 )^/2_ a .^}] 2 + e 2) 2 ' 

In the decoupled phase ( J r < J rc ) the impurity susceptibility diverges in the zero temperature limit (following the 
Curie law) while it vanishes in the screened phase. Since for rj^ = —0.2 the hybridization amplitude is smallest, 
the impurity-spin susceptibility becomes largest. Approaching the impurity QCP (J r — > J rc ), it shows a power-law 
divergence with an anomalous critical exponent of the ASL bulk. As shown in Fig. 3, such curves are well fitted with 

\imp ^ {Jr Jrc) ^ (20) 

where the scaling function is g(ri^) w 2 — 77^. It is valuable to consider how the behavior of the impurity susceptibility 
differs from that of the Pseudogap Kondo model[33[ which corresponds to the case of 77^, = 0. 

Next, we evaluate the impurity specific heat. The zero temperature formulation [Eq. (14)] for the impurity free 
energy can be transformed to the finite temperature version through the Wick rotation. Following Refs. [13|, [34[ , we 
find the impurity free energy at finite temperatures, 

f imp = N a r ^ 1—0(0 + 2 -^\b r ? - e rQx 

J -00 ePK + 1 Jr 



8 




FIG. 3: (Color online) Impurity spin susceptibility (log-log plot) as a function of the renormalized Kondo coupling for various 
ASL exponents. 




FIG. 4: (Color online) Specific heat coefficient (log-log plot) as a function of the renormalized Kondo coupling for various ASL 
exponents. 



with a rescaled temperature (3 r 1 = T r = T/A, where the "angle" function 0(f) is given by 

2 AT. 

+ 2 



8(?) ° -"( t [i + - + J + ? " - sig ° (0) torKI<1> 



^|!, r | J 8 in(^)«{« 2 -ir./2_^ r ,} 



tan- 1 ! — w ' ~T^FrT~T~7~~ : — ; — ) + -(1 -sign(£)) for |£ I > 1. (21) 



7T 



Here, the denominator and numerator in the angle function correspond to the real and imaginary parts of the kernel 
for the impurity free energy [Eq. (14)], respectively. 
We find the impurity entropy 



q dfi m p 



and specific heat coefficient 



C imv _ dS mp _N„ f 00 di e2 d ( 1 ^ 96(f) 



_ ^tmp _ ""imp _ -a I Z2.fi* _ \ )2l fnq\ 



Taking the zero temperature limit, we find the self-consistent results in Fig. 4, using the solutions of Eq. (15). The 
latter terms in Eq. (21) ensure that the impurity entropy is Si rnp /N a = In 2 in the decoupled phase, consistent with 
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FIG. 5: (Color online) Wilson's ratio (y axis-log plot) as a function of the renormalized Kondo coupling for various ASL 
exponents. 



our expectation. In the Kondo phase the impurity entropy becomes vanished. But, we note that more elaborate 
calculations result in small nonzero entropy contributions in the Kondo phase. [3ol [3a] The 7i mp coefficient shows a 
behavior similar to the impurity susceptibility Ximp, diverging as J r — > J rc . It exhibits the scaling behavior , 

with h(rj^) rs 2 — 0.2^.0 in our numerical analysis. 

Using the impurity susceptibility and specific heat coefficient, one can find the Wilson's ratio in the zero temperature 
limit 

W imp (T r - 0) = ^ . (25) 

Ximp T r >0 

In Fig. 5 we plot this value as a function of the rescaled Kondo coupling J r — J rc . Remember Wi mp — 2 in the 
Kondo effect of the Fermi liquid. Here, we also obtain a similar value for the Pseudogap Kondo model (rj^ = 0.001). 
An important observation is that the Wilson's ratio is strongly dependent on the ASL exponent. For the negative 
exponent the Wilson's ratio becomes enhanced while it gets suppressed for the positive one. This implies that the 
Wilson's ratio can be utilized as a probe for revealing the nature of SLs, more generally, criticality of the bulk system. 
This is the main message of the present study. 



IV. SUMMARY AND DISCUSSION 



It is valuable to remind several assumptions for solving the ASL Kondo problem. First of all, we have considered 
effects of staggered spin fluctuations on dynamics of a magnetic impurity, ignoring those of uniform spin correlations, 
since antiferromagnetic spin fluctuations are most susceptible in the ASL bulk, thus expected to give dominant con- 
tributions on this problem. In addition, ferromagnetic spin correlations do not show anomalous scaling, implying that 
such contributions would coincide with the Pseudogap Kondo effect, thus not so interesting. Our second assumption 
is in writing a spinon Green's function, where effects of gauge fluctuations are introduced in an anomalous critical 
exponent. In the present paper we have used the scaling exponent as a phenomenological parameter. Both assump- 
tions are compatible since the critical exponent of the staggered spin-spin correlation function is consistent with that 
of the single particle propagator when both critical exponents are evaluated in a gauge invariant way. 

The third one is rather an approximation than an assumption for solving the effective impurity action while the 
above two are basic assumptions for deriving the impurity action. In the slave-boson representation of this effective 
impurity action we have performed the large- N„ analysis introducing the hybridization order parameter. Although 
well-defined quasiparticle excitations do not exist in the case of ^ 0, it was shown that the I-QPT occurs between 
the local-moment state and the Kondo-screened phase. Evaluating the impurity spin susceptibility and specific heat 
coefficient, we found that the Wilson's ratio depends strongly on the ASL exponent. This result has an important 
physical meaning because the Wilson's ratio for the magnetic impurity reflects criticality of the bulk system. This 
conclusion will be available to general critical systems with exact screening particulary, where the expression of the 
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Kondo vertex is the same as that of the present paper. [37f In this respect the Wilson's ratio for the magnetic impurity 
may be one possible probe for measuring bulk criticality. 

It is interesting to compare the ASL Kondo physics with the Kondo effect in the Luttinger liquid, [38, H^] since the 
ASL can be considered as the high dimensional realization of the Luttinger liquid. In the Luttinger liquid the Kondo 
interaction term can be decomposed to the forward and backward scattering channels, analogous to the uniform and 
staggered ones in the ASL. It was shown that the forward scattering channel is irrelevant in the renormalization group 
analysis up to two- loop order. [3^| Similarly, in this paper we take into account only the antiferromagnetic correlation 
channel for the ASL Kondo effect, although it is not proven that the ferromagnetic channel is irrelevant. The backward 
impurity scattering in the Luttinger liquid was shown to cause anomalous scaling, in particular, power-law behavior 
of the Kondo temperature owing to the presence of the anomalous critical exponent in the Luttinger liquid. [39] This 
is basically the same as the ASL-Kondo effect that the ASL criticality results in anomalous scaling on the impurity 
physics, although there is no phase transition in the Luttinger liquid owing to one dimensionality. 

It should be noted that the present mean-field analysis is difficult to describe correct scaling behaviors in the 
over-screened phase since the hybridization order parameters are not regarded as dynamic variables but static ones. 
This approximation scheme seems to be more appropriate when quasiparticle excitations are well defined, thus the 
conventional Kondo screening picture is applicable. The slave-boson mean-field scheme can be improved using the non- 
crossing approximation, [35( where such hybridization parameters are taken to be dynamic variables, thus quantum 
fluctuations are more involved. Performing the Hubbard-Stratonovich transformation for the nonlocal (for time) 
hopping term in Eq. (9), we find an effective impurity action 

N„ , T AT,. 



s nca = I ^ xUiko -K + e x ) Xa + ^j- ^ K f 7ob s n - t x Q y 

a—I n—1 



dko 



where S x (fco) and Sh(fco) are fermion and boson self-energies, respectively, determined by the following self-consistent 
NCA-type equations 



(7=1 

E b (r - T >) = -T(t - r')(J2 ^M^(r')) 



with J-(k ) = ikf)F(kft). This kind of approximation is well known to catch non-Fermi liquid physics in the multi- 
channel Kondo model. [351 ] Scaling behaviors of both bosonic and fermionic self-energies are expected in the low energy 
limit, causing anomalous critical physics to this system even in the case of r\^, = 0. Inserting expected scaling forms for 
both the self-energies and renormalized Green's functions to the NCA equations, we would obtain the total anomalous 
scaling exponents which are expected to be sum of the critical exponents of the multichannel Pseudogap Kondo model 
and the ASL scaling dimension approximately, considering the presence of the ASL scaling exponent in T{t — r'). It 
will be interesting to examine how the scaling exponents in the conventional bulk are affected by the presence of the 
ASL exponent. 

Applying magnetic fields to the ASL, the impurity QPT is expected to disappear. Because external magnetic fields 
would result in finite density of states at the Fermi energy, the conventional Kondo physics may appear, where only 
the over-screened Kondo phase would occur, independent of the Kondo interaction. Furthermore, gauge fluctuations 
would be dissipative due to the finite density of states, and the bulk system becomes more "Fermi liquid" -like, 
supporting the above expectation. 

In the present analysis we did not consider scattering due to randomly distributed disorder potentials. One of 
the present authors has studied the role of random potentials in the ASL, and found that such a spin liquid phase 
remains stable against weak disorders because massless Dirac spinons at the interacting fixed point live in higher 
spatial dimensions than two owing to the presence of the anomalous critical exponent. [401] Remember the presence 
of the delocalization transition above two spatial dimensions. However, it is not clear whether the diffusive nature 
appears or not in the ASL. If so, the presence of finite density of states due to random potential scattering may 
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destroy the I-QPT as the case of magnetic fields. 



[1] R. Coldea, D. A. Tennant, A. M. Tsvelik, and Z. Tylczynski, Phys. Rev. Lett. 86, 1335 (2001); R. Coldea, D. A. Tennant, 

and Z. Tylczynski, Phys. Rev. B 68, 134424 (2003). 
[2] Y. Shimizu, K. Miyagawa, K. Kanoda, M. Maesato, and G. Saito Phys. Rev. Lett. 91, 107001 (2003); Y. Kurosaki, Y. 

Shimizu, K. Miyagawa, K. Kanoda, and G. Saito Phys. Rev. Lett. 95, 177001 (2005); A. Kawamoto, Y. Honma, and K. I. 

Kumagai, Phys. Rev. B 70, 060510 (2004). 
[3] J. S. Helton, K. Matan, M. P. Shores, E. A. Nytko, B. M. Bartlett, Y. Yoshida, Y. Takano, A. Suslov, Y. Qiu, J.-H. 

Chung, D. G. Nocera, and Y. S. Lee, Phys. Rev. Lett. 98, 107204 (2007); P. Mendels, F. Bert, M. A. de Vries, A. Olariu, 

A. Harrison, F. Due, J. C. Trombe, J. S. Lord, A. Amato, and C. Baines, Phys. Rev. Lett. 98, 077204 (2007); Oren 

Ofer, Amit Keren, Emily A. Nytko, Matthew P. Shores, Bart M. Bartlett, Daniel G. Nocera, Chris Baines, Alex Amato, 

|cond-mat/0610540"| 

[4] Marc Bocquet, Fabian H. L. Essler, Alexei M. Tsvelik, and Alexander O. Gogolin, Phys. Rev. B 64, 094425 (2001). 
[5] S. V. Isakov, T. Senthil, and Y. B. Kim, Phys. R ev. B 72, 174417 (2005). 
[6] Yi Zhou and Xiao- Gang Wen jcond-mat /0210662~j 

[7] J. Alicea, O. I. Motrunich, and M. P. A. Fisher, Phys. Rev. Lett. 95, 247203 (2005); J. Alicea, O. I. Motrunich, and M. P. 

A. Fisher, Phys. Rev. B 73, 174430 (2006). 
[8] T. Senthil's talk in ICTP. 

[9] S. Florens, L. Fritz, and M. Vojta, Phys. Rev. Lett. 96, 036601 (2006). 
[10] M. Vojta, C. Buragohain, and S. Sachdev, Phys. Rev. B 61, 15 152 (2000). 

[11] Ying Ran, Michael Hermele, Patrick A. Lee, and Xiao-Gang Wen, Phys. Rev. Lett. 98, 117205 (2007). 
[12] J. B. Marston and I. Affleck, Phys. Rev. B 39, 11538 (1989). 

[13] D. Withoff and E. Fradkin, Phys. Rev. Lett. 64, 1835 (1990); C. R. Cassanello and E. Fradkin, Phys. Rev. B 53, 15079 
(1996). 

[14] L. Fritz and M. Vojta, Phys. Rev. B 70, 214427 (2004); M. Vojta and L. Fritz, Phys. Rev. B 70, 094502 (2004). 
[15] K. Ingersent, Phys. Rev. B 54, 11936 (1996). 

[16] M. Hermele, T. Senthil, M. P. A. Fisher, P. A. Lee, N. Nagaosa, and X.-G. Wen, Phys. Rev. B 70, 214437 (2004). 

[17] M. B. Hastings, Phys. Rev. B 63, 014413 (2001). 

[18] D. H. Kim and P. A. Lee, Annals Phys. 272, 130 (1999). 

[19] T. Appelquist, D. Nash, and L. C. R. Wijewardhana, Phys. Rev. Lett. 60, 2575 (1988). 

[20] The renormalized spinon propagator is expressed as G^g L (k) = G c ^ 1 (fc) + E(fc), where G ( 7 1 (fe) = i7 M fc M is the inverse of the 
bare spinon propagator, and E(fc) the spinon self-energy resulting from long-range gauge interactions. In the 1/N a expansion 

the self-energy is represented as £(fc) = J ^iiGo{k + g)7„D M „(g), where D llv {q) is the renormalized propagator of the 

U(l) gauge field due to polarization of massless Dirac fermions. The gauge propagator is obtained to be D llv {q) « n~„ (q) = 

(d/iv ~ ^a" ^ in the Lorentz gauge, where Tl^ v (q) — N a J ^p-Tr[Go(fc)7 M Go(fc + q)jv] is the polarization function of 

Dirac fermions. Inserting this gauge propagator into the expression of the self-energy, one can find the spinon self-energy 

of logarithmic momentum dependence E(fc) = irj^j^k^lnyj^ , where 77,/, is the anomalous exponent and A, the momentum 

cutoff. The absolute value of the exponent rj^ is proportional to the inverse of the flavor number, i.e., |tj^| ~ N^ 1 . 
Such logarithmic momentum dependence should be considered as the lowest order in the algebraic function. As a result, 
one can obtain the algebraically decaying spinon propagator [Eq. (2)] from the following nonperturbative consideration 

GaslW = il»K [l + ^n( |]] »t7„fc„(£)"* . 

[21] W. Rantner and X.-G. Wen, |cond-mat/0105540| 

[22] J. Ye, Phys. Rev. B 67, 115104 (2003). 

[23] D. V. Khveshchenko, Phys. Rev. B 65, 235111 (2002). 

[24] V. P. Gusynin, D. V. Khveshchenko, and M. Reenders, Phys. Rev. B 67, 115201 (2003). 

[25] M. Franz, Z. Tesanovic, and O. Vafek, Phys. Rev. B 66, 054535 (2002). 

[26] P. M. Chaikin and T. C. Lubensky, Principles of condensed matter physics, (Cambridge University Press, Cambridge, 

England, 1995). 

[27] Y. Ran and X.-G. Wen, |arXiv:cond-mat/ 0609620 

[28] M. Hermele, T. Senthil, and M. P. A. Fisher, Phys. Rev. B 72, 104404 (2005). 

[29] A. Tanaka and X. Hu, Phys. Rev. Lett. 95, 036402 (2005). 

[30] M. Franz, T. Pereg-Barnea, D. E. Sheehy, and Z. Tesanovic, Phys. Rev. B 68, 024508 (2003). 

[31] W. Rantner and X.-G. Wen, Phys. Rev. B 66, 144501 (2002). 



12 



[32] The saddle-point equations of the Pseudogap Kondo model can be obtained from Eq. (14) in the limit of 77^ — » 

1 _ 1 r d ,f 2 Mi+^)[i+^i^i 2 Mi+^)] 



Jr 27T 2 

Qx = _l 

No- 7r./ n T 2ri , Jv^, 



dx- 



giving rise to Eq. (16) for the impurity QCP. 
[33] In the Pseudogap Kondo model the impurity susceptibility is given by 



n J ( X 2[ 1 + ^| 6r | 21n ( 1+ + 

Remember that the impurity susceptibility based on the large- N a analysis exhibits a logarithmic correction due to the 

upper critical dimensionality. 13] 
[34] S. Doniach and E. H. Sondheimer, Green's Functions for Solid State Physicists, (The Benjamin/Commings Publishing 

Company, Massachusetts, USA, 1982). 
[35] D. L. Cox and A. E. Ruckenstein, Phys. Rev. Lett. 71, 1613 (1993); O. Parcollet, A. Georges, G. Kotliar, and A. Sengupta, 

Phys. Rev. B 58, 3794 (1998). 
[36] O. Parcollet and A. Georges, Phys. Rev. Lett. 79, 4665 (1997). 

[37] Note that the 0(3) critical bulk described by spin 1 critical fluctuations is not described by the present Kondo model. 
The present results cannot be applied to all critical bulk systems in spite of their generality because the vertex in Kondo 
interaction differs from each system, depending on quantum numbers of critical modes in the bulk system. 

[38] D.-H. Lee and J. Toner, Phys. Rev. Lett. 69, 3378 (1992); A. Furusaki and N. Nagaosa, ibid. 72, 892 (1994). 

[39] M. D. Kim, C. K. Kim, K. Nahm, and C.-M. Ryu, J. Phys.: Condens. Matter 13, 3271 (2001). 

[40] Ki-Seok Kim, Phys. Rev. B (R) 70, 140405 (2004); Ki-Seok Kim, Phys. Rev. B 72, 014406 (2005). 



